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Abstract 

We perform a comprehensive analysis of the dynamical magnetic susceptibility xCQi ^) the 
slave-boson mean-field scheme of the bilayer t-J model. We use model parameters appropriate 
for YBa2Cu306+a; (YBCO), a typical bilayer high- Tc cuprate compound well studied by neutron 
scattering experiments. In the d-wave pairing state, the strongest magnetic spectral weight appears 
at q = Q = (vr, vr) and u = uJq^, and spreads into a diamond-shaped shell around Q in q 
space for lu < uJq^. This weight is due to a collective mode, namely a particle-hole bound state, 
which has a downward u versus q dispersion around Q. Within the high intensity shell, the 
incommensurate (IC) signals at q = (tt, vr it Inr]) and (vr it 2vrr/, vr) tend to be stronger than 
the diagonal incommensurate (DIG) signals at q = (vr it 2vrr/', vr it 27rrj'), especially for a large 
hole density 5. For to <^ uJq^ the IC signals completely disappear and the weight remains only 
around the DIC positions. For lu > i^q*^ strong signals of Imx(q, uj) tracing an upward dispersion 
are found and interpreted as an overdamped collective mode near Wq^. In the normal state, 
Imx(q, w) has a broad peak at q = Q. That is, the IC and DIC signals appear only in the d- 
wave pairing state. We also study effects of a small orthorhombic anisotropy, which is intrinsic in 
untwinned YBCO crystals. Because of electron-electron correlations favoring d-wave shaped Fermi 
surface deformations (dFSD), we expect an enhanced anisotropy of magnetic excitation spectra. 
This effect is particularly pronounced for low 5 and at relatively high temperature. The present 
theory provides a rather detailed microscopic explanation of the most salient properties of magnetic 
excitations observed in YBCO. 

PACS numbers: 74.25.Ha, 74.72.Bk, 74.20.Mn, 71.10.Fd 
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I. INTRODUCTION 



The undoped parent compounds of high-Tc cuprate superconductors are antiferromag- 
netic Mott insulators. With carrier doping the antiferromagnetism disappears and a super- 
conducting state with a high transition temperature is reahzed. Antiferromagnetic (AF) 
correlations however remain even in the superconducting state and the understanding of 
magnetic properties has been widely recognized as a major issue in the theory of high-Tc 
cuprates. 

Magnetic correlations are directly measured by inelastic neutron scattering experiments. 
Most of the experiments were performed for the cuprate superconductors YBa2Cu30y 
(YBCO) and La2_xSr^.Cu04 (LSCO). Early studies suggested that YBCO had commen- 
surate antiferromagnetic correlations, that is, the peak signal in the imaginary part of 
the dynamical magnetic susceptibility xi^j ^) various frequencies appeared at q = 
Q = (tt, 7r),[l| while LSCO had incommensurate (IC) correlations in the sense that the 
peak shifted to q = (tt ± 2717?, vr) and (vr, tt ± 2777?), where ri parametrizes the degree of 

n n n n 

incommensurability. 12, |3|, |^ |5|{ In 1998, measurements with improved resolution revealed 
that the seemingly commensurate signal of YBCO was actually composed of four peaks at in- 
commensurate wavevectors. ol Indications for an incommensurate signal in YBCO appeared 
already in several earlier observations, 13, M but remained ambiguous due to a limited resolu- 

nn 

tion. The peak positions were identified as q = (7r±27r?7, vr) and (vr, 7r±27rr7),[9|,|10[ the same 
geometry as in LSCO. The IC peaks, however, appeared only below the superconduct ing 
transition temperature Tc or possibly below the pseudogap temperature T*, 

mum 

in sharp contrast to LSCO, where the IC signals were observed at temperatures well above 

The size of the incommensurability rj depends on the excitation energy u. For YBCO, r) 
decreases with increasing u and vanishes at a specific frequency uJq^. The commensurate, 
so-called "resonance peak" at w^was extensively investigated for optimally doped 13 
0, Q and underdoped samplesQ, Q, Q, Q before the observation of the IC signals 
in YBCO. The resonance peak is now regarded to be continuously connected with the 
IC signals observed at lower frequencies, that is the peak disperses smoothly downwards 
to lower frequencies when q is shifted away from Q. jiol. 111 . 22 1 Above the resonance 



energy, possible IC structures were observed. JjQ] It was recently found that the strongest 
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weight appeared at diagonal incommensurate (DIG) positions q = (vr ± 2711], vr ± 2711]). 23| 
On the other hand, for LSCO the energy dependence of rj was relatively weak compared 
to YBCoQ and such a robustness of IC structures was often contrasted with the behavior 
in YBCO. However, recent high energy neutron scattering data for La2_2,.Ba2^Cu04 with 



X = 0.125 revealed t 
peak around 55meV. 



lat the IC peaks dispersed with u and merged into a commensurate 



24| Above this energy the data showed weak DIG peaks as in YBGO. 



Early theoretical work[ 



Q y, 27, Q pointed out that the differences of magnetic ex- 
citations between YBGO and LSGO could be traced back to the difference of the Fermi 
surface (FS) shape, which we refer to as the fermiologv scenario. The FS difference was 



indeed predicted by LDA band calculations 29|, |30|, |3l| and angle-resolved photoemission 



spectroscopy (ARPES) supported the LDA prediction for YBGO;j3al33l ARPES data for 
LSGO were not available for a long time. After the experimental observation of IG peaks in 



YBGO, further detailed calculations 



ionsjsj. 
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37 



38, 



3S 



within the same theoretical 



framework as in the early days showed that essential features of magnetic excitations in 
YBGO were captured in the fermiology scenario. 

Fermiology theories for LSGO systems were running into problems when ARPES data. 



first reported in 1999, 



4l| suggested that the FS of LSGO looks similar to that of YBGO in 



the doping region below x ~ 0.15 .fA^ If so, fermiology scenarios predict that magnetic exci- 
tations should be essentially the same in both YBGO and LSGO, in contradiction with the 
experimental data. This problem was considered within the slave-boson mean-field scheme 

Qn 
It was shown [44] that the model has tendencies toward orientational 

symmetry breaking of the square lattice symmetry, leading to a d-wave shaped FS deforma- 
tion ((iFSD): the FS expands along the axis and shrinks along the ky axis or vice versa. 
Assuming a coupling to the low temperature tetragonal lattice distortions or their fluctu- 
ations, one could understand both the observed FS shape and the structure of magnetic 
excitations in LSGO consistently. 3, 42| 

An essentially different picture for magnetic excitations was proposed in 1995, j3| the 
spin-charge stripe scenario, according to which IG magnetic excitations are mainly con- 
trolled by charge-stripe order tendencies in the Gu02 plane. This scenario is based on 
the observation of a charge order signal whose wavevector is just twice as large as the IG 
magnetic wavevector. Although charge order was observed only in a few high-Tc cuprate 
materials with specific doping rates, |3, li^ 1^ and in addition the signal was rather 
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weak, the spin-charge stripe scenario attracted much interest. [5^ 

The spin-charge stripe scenario predicts one-dimensional magnetic signals, that is, IC 
magnetic peaks should appear either at q = (tt ± 27r?7, vr) or (tt, tt ± 2TTr]). Such a one- 
dimensional pattern was actually inferred from experiments for partially untwinned YBCO 
crystals 53| as a strong support for the stripes. However, subsequent neutron scattering 
studies for almost fully untwinned YBCOjSJ, |55| revealed that magnetic excitations were 
two-dimensional and had four IC signals. The IC pattern however lost fourfold symmetry 
around q = Q and exhibited some anisotropy between the and qy direction. There are 
several possible explanations for this anisotropy. (i) The bare band structure effect, due 



to the orthorhombicity of the crystal structure 



56|, |57| or the CuO chains in YBCO.|58| 



(ii) The tendency to a rfFSD due to electron-electron correlations, which enhances the bare 
a,.s„t.„py.Q3 (Hi) Effec. of cKa.,e st.pe TKe putative charge ..pes 

should be strongly fluctuating to be consistent with the experimental observation that the 
magnetic signals form a two-dimensional anisotropic geometry. 

We focus on the possibility of an enhanced anisotropy due to correlations favoring a rfFSD. 
The (iFSD tendency is generated by forward scattering interactions of electrons close to the 
FS near (vr, 0) and (0, vr). It was first found in the t-J model j^] and Hubbard model, 



and then analyzed in more detail 62 



subsequent works. 



68, 



6S 



TOl, EJlf the (iFSD tendency is very strong, it can lead to a 



3, 



64, 



66 



67l | and also for other models in various 



spontaneous breaking of the orientational symmetry of the Fermi surface from tetragonal to 
orthorhombic. Referring to a stability criterion for Fermi liquids by Pomeranchuk. j72l| some 
authors termed this instability "Pomeranchuk instability" . 



The (iFSD competes with (i-wave singlet pairing 
stripes, j4 



. H H 0, E , 65 , 3 It differs from other 



frequently discussed ordering tendencies such as antiferromagnetism, spin-charge 



ItsI i staggered flux. j76j| and (i-density wave, [77] which are driven by interactions 
with a large momentum transfer near q = Q, while the (iFSD is generated from forward 
scattering. The (iFSD breaks the orientational symmetry of the square lattice and has the 
same reduced symmetry as the electronic nematic phase proposed by Kivelson, Fradkin and 
Emery. [3] Their route to this phase is however not that of a Fermi surface instability, but 
rather via partial melting of stripe order. 

To discuss high- Tc cuprates we should treat both the (iFSD and the (i-wave singlet pairing 
on an equal footing. We use the slave-boson mean- field scheme of the t-J model; both 



4 



tendencies are generated by the J term. The dFSD tendencies are suppressed by the (i-wave 
singlet pairing instabihty such that no spontaneous Fermi surface symmetry breaking takes 
place. 4JI However, significant (iPSD correlations survive. 4J, |66[ This correlation effect can 
drive a sizable enhancement of the Fermi surface anisotropy by a small bare band anisotropy. 
We analyze this dFSD effect on magnetic excitations by using the anisotropic t-J model. The 
idea that dFSD correlations may enhance the in-plane anisotropy of magnetic excitations in 
YBCO has been pursued already in a recent work by Kao and Kee 3] ■ In that work ciFSD 
correlations and pairing are taken into account via a suitable ansatz for the renormalized 
band structure and the superconducting gap function, while we actually compute these 
interaction effects, starting from a microscopic model. 

Many features of magnetic excitations in YBCO are captured already by the (usual) 



isotropic t-J model 



model 
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Vlagnetic excitations in YBCO are well characterized by the bilayer 
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80 



and therefore have odd and even channels. While most 
neutron scattering studies were confined to the odd channel so far, recent experiments suc- 
cessfully detected the even channel also. We perform a comprehen- 
sive analysis of magnetic excitations in the bilayer t-J model including the above mentioned 
(iFSD effects. We show that prominent features of magnetic excitations of YBCO are well 
captured in this framework and confirm several theoretical insights obtained earlier. The 
(iFSD effect turns out to provide a natural scenario to understand the observed anisotropy 
of magnetic excitations, especially at low doping and relatively high temperature. Together 



with previous studies of ma gne tic excitations in the single layer t-J model 



I45I and 



in the bilayer t-J model, \ the present study provides a comprehensive understanding 
of magnetic excitations in YBCO from the fermiology viewpoint. Since the antibonding FS 
can be easily deformed to be open by c/FSD correlations, our work is also relevant for the 
interpretation of ARPES data for untwinned YBCO. 

The article is structured as follows. In Sec. II, we introduce the bilayer t-J model and 
present the slave-boson mean-field scheme. In Sec. Ill, we discuss self-consistent mean-field 
solutions of the bilayer system. Sec. IV is dedicated to magnetic excitations. After a brief 
qualitative discussion of the dynamical magnetic susceptibility, we present numerical results 
in two parts: (i) the isotropic case, where comprehensive results are provided to discuss 
prominent features of magnetic excitations in YBCO, and (ii) the anisotropic case, where 
the (iFSD effect on magnetic excitations is investigated. In Sec. V, we summarize our results 
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through a comparison with experimental data, and conclude in Sec. VI. 

II. MODEL AND FORMALISM 

We analyze the bilayer t-J model on a square lattice 

H = — ^ trr'cl^Cr'a + ^ Jrr' Sr " Sr' (l) 

r,r',(T (r,r') 

defined in the Fock space with no doubly occupied sites. The operator cj^ (cro-) creates 
(annihilates) an electron with spin a on site r, while Sr is the spin operator. The site 
variable r = (x, y, z) runs over the bilayer coordinates, that is (x, y) indicates a site on 
the square lattice and each layer is denoted by z = or 1. We assume periodic boundary 
conditions so that A;^ = or tt. Jrr' (> 0) is a superexchange coupling between the nearest 
neighbor sites along each direction, x, y, and z. We take into account hopping amplitudes 
trr' between r and r' up to third- nearest neighbors. We denote Jrr' and trr' by using the 
conventional notation defined in Fig. ^ 

We introduce the slave particles, f^a and for, as Cro- = h\.fra, where f^a (K) is a fermion 
(boson) operator that carries spin a (charge e), and Sr = ^fra^a/sfr^ with the Pauli ma- 
trices cr = [(7^,(7^,(7^). The slave bosons and fermions are linked by the local constraint 
blK + Yla frafra = 1- Thls Is an cxact transformation known as the slave-boson formalism. 
We then decouple the interaction with the so-called resonating-valence-bond (RVB) mean 
fields:0 XT={Y.afLfr'a), {blK'), and A^={fr^fr'i - /rj/r't), with r = r' - r. These 
mean fields are assumed to be real constants independent of sites r. We approximate the 
bosons to condense at the bottom of the band, which is reasonable at low temperature and 
leads to (blK') = 6, where 6 is the hole density. The resulting Hamiltonian reads 

with a global constraint J^aifrafra) = ^ — the k summation is over \k^(^y-j\ < n, and 
kz = and vr. Here = eli + ej^ — f^, with the in-plane (c-axis) dispersion ef, (ej^), and 
Ak = a||. + Ak"*", with the singlet pairing gap in (out of) the plane A||, (Aj^); fi is the chemical 
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potential. The explicit momentum dependence of the dispersion is given by 
4 = -2 (^J + ^JxXx^ cos kx + (^y6 + ^JyXy] cos ky 

+2t'6 cos k^ cos ky + cos 2kx + ty6 cos 2ky 
3 \ 

^±0^ + -^JzXz I + 2 (t_La;5 COS fca;. + t^y5 COS 



+ 2i'_i_(5 COS kx COS /cy + t"^^5 cos 2 /c^; + ^5 cos 2ky) 



cos kz , 



and that of the gap function by 



(3) 



(4) 



(5) 
(6) 



Al = -- {JxA^ cos k^ + JyAy COS ky) , 
3 

Aj^ = --JzAz coskz . 



In the bilayer high-Tc cuprates, a (cos kx—cos k^-type c-axis dispersion was computed 
and actually observed in Bi2Sr2CaCu2O8+5.0, Hi To adapt Eq. (jH) to this behavior, we 
choose transverse hopping amplitudes as follows: t±o = '^^^^^tz-, t±x = t±y = 0,t'j_ - 

2 2 

^±x ~ and t'^y = -^tz- The parameters 'jx and 7^ allow for a convenient parametrization 
of an in-plane anisotropy as specified below. The resulting c-axis dispersion reads 



2 ' 



2 3 

tzS {-fx COS kx - Jy COS ky) +-JzXz 



COS k. 



(7) 



which has the expected form in the isotropic case 73^. = 7^ = 1. 

YBa2Cu30j^ has an orthorhombic crystal structure for y > 6.4, where the supercon- 
ducting state is realized at low T. 9^ Such orthorhombicity yields xy-anisotropy, which we 
incorporate by introducing a single parameter a as follows: 

tx=t{l + a/2), ty = t{l-a/2), 
Jx = J{l + a), Jy = J{l-a), 



t:=t"{l + a/2), t;' = t"(l-a/2), 
7a; = l + a/4, 7j^ = l-a/4. 



(8) 
(9) 
(10) 

(11) 



With this parametrization, tx{y), i'x{y)i and i"j^x(y) have the same degree of x?/-anisotropy, 
while the anisotropy of Jx{y) is twice as large, as imposed by the superexchange mechanism. 
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We determine Xt and with r = x, 2; by solving the following self-consistency equa- 
tions numerically: 




where = a/^^ + ^ total number of (bilayer) lattice sites. 



III. SELF-CONSISTENT SOLUTION 



The material dependence of hig h- T^ cu prates is mainly taken into account by different 
choices of band parameters. j3, 0,0,13 We use the following parameters for YBCO: 

t/J = 2.5, t'/t = -0.3, t"/t = 0.15, J,/ J = 0.1, tjt = 0.15 , and a = -0.05 . 



This choice has been done judiciously. Ab initio calculations [9^ indicate that realistic values 
for t/J lie in the range 2 — 5. Even within this restricted interval, we found that the behavior 
of Imx(q, uj) strongly depends on the choice of t/ J. We have chosen t/J in such a way that 
the energy of the resonance mode at q = (tt, vr) is maximal at optimal doping, that is near 
0.15 (see Fig. [rUj) . to roughly agree with experiments 
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at optii 
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8J The ratios t'/t 
8_7ll The value of Jz/ J is 
9q The bilayer system 



and t" /t are extracted from an LDA band calculation for YBCO. 
fixed rather uniquely by the optical magnon energy in YBCO. [95, 
has two FSs, the bonding band FS and the antibonding band FS. With increasing 6, the 
latter can change from a hole-like FS to an electron-like FS in the region we are interested 
in. The hole density at which such a FS topology change occurs strongly depends on tz- 
With our choice of t^/t, which is a bit larger than LDA estimates, 83] the antibonding FS 
becomes electron-like for 6 > 0.20. Some ARPES data for Bi^Sr?CaCu?iOs+,s indeed suggest 



89 



ARPES data for 



an electron-like antibonding FS for sufficiently high overdoping. 
YBCO with such a high doping are not yet available. 

It is not easy to determine the anisotropy parameter a; even its sign is not obvious. 
The orthorhombic YBCO has a lattice constant anisotropy a < b, where a{b) is the lattice 
constant along the x{y) direction. This contributes to an enhancement of t^- On the other 



hand, YBCO contains CuO chains along the y direction. A small hybridization with the 
chain band enhances ty, an effect opposite to the lattice constant effect. If one estimates 
a band anisotropy from the predicted FSs in Ref. js^], a small negative a is obtained. A 
most recent LDA calculation supports this estimate and predicts that ty is larger than 
by about 3 — 4%.j9^ Hence we assume < ty and set a = —0.05 when analyzing the 
in-plane anisotropy of Imx(q, u). We will also present results for the isotropic case, a = 0, 
for comparison. 

We introduce the following convenient notation: 

Xx ~l" Xy Xx Xy fi c\ 

Xo = Xd = (15) 

A, = ^^, A,= \ (16) 

where Xd is an order parameter of the (iFSD, and Ad{As) is the d-wave (s-wave) in-plane 
singlet pairing amplitude; we choose A^,. > 0. 
For a = 0, there is no (iFSD, 



4i|66| d-wave pairing state 

is stabilized at low T, that is, A^ = A^ = 0. Note that a finite A^ contributes to an s-wave 
pairing component of A^ [see Eq. ©]. The temperature dependences of Xo? Xz, and A^ are 
shown in Fig.|21^a) and their 6 dependences in Fig. |2fb). Xo ^-^id Xz depend very weakly on T. 
In particular, the onset of does not affect xo and Xz substantially. At 5 = 0, Xz becomes 
zero, that is, the layers are decoupled regardless of the interlayer coupling in the present 
mean-field framework, where antiferromagnetic order is not taken into account. The layers 
couple weakly only for finite 6. The d-wave pairing amplitude Ad increases with decreasing 
6, which however should not be interpreted as an enhancement of a superconducting order 
parameter at lower 6, but rather as the increase of a one-particle pseudogap energy scale in 
the slave-boson mean-field scheme. 0| 

When the anisotropy parameter a is introduced, the system loses tetragonal symmetry, 
and Xd, As, and A^ can become finite. Since the temperature and doping dependences of 
Xo, Ad, and Xz for a = —0.05 are almost the same as in Fig. HI we focus on Xd, Ag, and 
Az- The temperature dependences of Xd, Ag and A^ are shown in Fig. |2fa). Xd increases 
with decreasing T, which is due to the development of dFSD correlations, and exhibits a 
cusp at the onset of singlet pairing, which is denoted as Trvb in the slave-boson theory; 
note that A^^, A^,, and Ag have the same onset temperature. Below Trvb, (iFSD correlations 
are suppressed leading to a suppression of Xd- But Xd saturates at a finite value at zero 



temperature and is still enhanced compared to Xd at high T. Since the renormalized band 
is mainly characterized by the renormalized nearest-neighbor hopping ir = tr5 + '^JtXt 
[Eq. (jni)], this enhancement contributes to the band anisotropy. Defining At = \tx — iy\/2 
and to = iix + ^j/)/2, we plot 2At/to in Fig. IHl^b) for several choices of 5. Although the bare 
anisotropy is only 5%, the renormalized band anisotropy is enhanced up to around 25% for 
5 = 0.08 by (iFSD correlations. The enhancement is less prominent at higher 6. 

Figure Efc) shows the doping dependences of and A^ at low T. Although A^ increases 
at lower 6, its magnitude remains very small compared to A^. Note that we obtain a positive 
As for a = —0.05, that is A^. > \Ay\ for < Jy, which means that the s-wave component 
reduces the anisotropy of |A||.| caused by the anisotropy of J [see Eq. The out-of-plane 
pairing A^ becomes finite at finite 6. Its sign is determined uniquely through the linear 
coupling to Ad- It becomes positive for the present band parameters. In the antibonding 
band {k^ = vr), A^ enhances the magnitude of A^ at (vr, 0) and reduces it at (0, vr), and vice 
versa in the bonding band {kz = 0). 

The FS at low T, which we define by = in the singlet pairing state, is shown in 
Figs. El^a)-(c) for several choices of 6. The bonding band {kz = 0) forms the outer FS and 
the antibonding band {kz = vr) the inner FS. For a = 0, the inner FS changes into an 
electron-like FS for 6 > 0.20 while the outer FS stays hole-like. For a = —0.05, the inner FS 
becomes open for all 6 in Fig. EJ The inner FS always opens already for a tiny anisotropy 
parameter for doping around 6 ~ 0.20, since the original inner FS is close to (tt, 0) and 
(0, vr) in this case. It depends sensitively on the size of a whether an open FS is realized at 
low 6 where ciFSD correlations become large but the original inner FS is away from (tt, 0) 
and (0, tt). The antibonding band and the bonding band have different hole densities, 6a 
and 6b, respectively, as shown in Fig. From the inset of Fig. El we see that the antibonding 
band has a hole density more than 1.5 times higher than the total (average) hole density 6. 

IV. MAGNETIC EXCITATIONS 

Next we investigate magnetic excitation spectra in the bilayer t-J model using the self- 
consistent mean-field solutions. The irreducible dynamical magnetic susceptibility Xo(ci; 
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is given by 

xo(q, 



1 y ( tanh ^ - tanh ^] 



+ 2 k+q 



- -Ek+q + + ir 



tanh ^ + tanh — 



2T 



2T 



X 



+ 



1 



. Eli + -£^k+q + + ir -Ek + -E'k+q 

where F is a positive infinitesimal and 



ir 



(17) 



^k,k+q 



1 A ^ a^k+q + AkAk+q A 



(18) 



-f'k-f'k+q 

Note that the fcz-component of k is summed only over the two values kz = and it, cor- 
responding to a bonding and an antibonding band, respectively. Particle-hole scattering 
processes are therefore intraband processes for Qz = and interband processes for = vr. 
The former is called the even channel and the latter the odd channel. 

In a renormalized random phase approximation (RPA) 2^, the dynamical magnetic 
susceptibility x(ct) ^) is given by 



x(q, ^) 



xoiq, 



where 



J(q) = 2r( J^. cos + Jy cos qy) + Jz cos 



(19) 



(20) 



with a renormalization factor r. In the plain RPA one has r = 1, which leads to a divergence 
of x(q) 0) at q ~ (tt, tt, tt) in a wide doping region (5 < 0.20), signaling an instability 
toward the AF state. However, several numerical studies of the t -J model indicate that 
the an«fe„o,„agne* i,..abn.y ,s ove..e.ti.a.ed by the RPA.Q U U Fluctuations 
not included in RPA obviously suppress the instability. This can be taken into account 
in a rough and phenomenological way by setting r < 1 . 0, 0| Here we choose the value 
r = 0.5, which confines the AF instability to 5 < 0.064, consistent with actual YBCO 



samples. |9lll02| 



In the following we will frequently specify by refering to the "even" or "odd" channel, 
and include only the remaining and g^ variables in q. 
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A. Basic properties 



Before calculating the dynamical magnetic susceptibility numerically, we first discuss 
generic properties of Imx(q, uj), which hold in both odd and even channels. At zero tem- 
perature the (i-wave superconducting state is realized and 



(21) 



Hence, Imxo(ci? ^) has a threshold energy defined by 

<:^th(q) = min{Ek + -Ek+q, k G BZ} 



(22) 



above which continuum excitations start. ci;th(q) is sketched by the solid line in Fig. Inland 
has a maximum at q = Q. Along the diagonal direction q = :^(5', q), gapless excitations 
appear, which are due to scattering processes between the (i-wave gap nodes. 

The denominator in Eq. ()19|) vanishes if both real and imaginary parts vanish, that is 



1 + J(q)Rexo(q, <1 = 0, 
Imxo(q, <') = 0. 



(23) 
(24) 



As widely discussed in the literature. 



34 



36 



38, 



3S 



40 



58, 



103 



104, 



105| Eqs. (El and 



can be simultaneously satisfied when Imxo(q; ^) jumps from zero to a finite value at 
Lj = u;th(q), which gives rise to a logarithmic divergence of Rexo(q) such that there is 



always a solution u 



< Uth{(i)- This solution describes an ingap collective mode of 



particle- hole excitations with spin 1 and charge zero. Note that q is not restricted to Q. In 
fact, the collective mode appears also around q = Q , in particular at IC or DIG wavevectors 
as sketched in Fig. El Expanding J(q)Rexo(q? ^) around q = Q and u = Uq', we obtain 
the asymptotic form of the dispersion 



UJ. 



Q 



+ 



where 



ra 



[jm 



.1 9V(q)Rexo(q, ^) / 5Rexo(q, ^) 



dui 



(25) 



(26) 



q=Q, 1^=0;^'' 



The denominator in Eq. ()26|1 is usually positive for uj < co'th(Q)- The numerator is expected 
to be negative close to the commensurate AF instability since | J(q)Rexo(q; uj)\ has a com- 
mensurate peak and J(q) < 0, leading to an upward dispersion (dotted line in Fig. P). 
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With hole doping, the peak of |J(q)Rexo(q, ^)\ usually shifts to IC positions, and thus its 
curvature at q = Q becomes positive, that is, c^q*^ has a downward dispersion (bold solid 
line in Fig. IHI). When the collective mode mixes with the continuum, it disappears due to 
overdamping. 

The infinitesimal F in Eq. (|17|) is replaced by a finite small value in numerical calculations, 
leading to a cutoff of the logarithmic divergence of Rexolq, A similar cutoff occurs also 
for finite T. Then Eq. fl23|) can be satisfied only when Imxolq, i^) has a sufficiently large 
jump at u = cc'th(q) or when Rexo{(i, i^) is sufficiently large; the latter is always the case 
near the AF instability. Practically therefore the collective mode is well defined only in a 
limited momentum region (Fig. IHI). Moreover, since Imxo(q; ^) becomes finite for a finite 
F or a finite T, the collective mode has a finite life time, which yields a finite width for the 
peak in Imx(q, u). 

B. Numerical calculation 

Now we perform extensive numerical calculations of Imx(q, u) with a small finite F = 

0. 01J in the denominator of Eq. p7|) . Although the choice of a finite F is mainly done 
for numerical convenience, it also simulates damping of electrons by static defects in real 
materials, and broadening due to limited energy resolution in inelastic neutron scattering 
experiments. 

This section is composed of six parts, where the first three deal with the isotropic case 
[a = 0) and the last three with the effects due to anisotropy {a 7^ 0). The former are 
intended to give a comprehensive analysis of magnetic excitations in bilayer cuprates such as 
YBCO, without taking effects which are specifically due to in-plane anisotropy into account: 

1. (q, uj) maps of Imx(q, uj) for a sequence of 6 for both odd and even channels, 2. q maps 
for a sequence of uj, and 3. T dependence. The latter are dedicated to anisotropy effects on 
magnetic excitations and their enhancement by rfFSD correlations: 4. q maps for a sequence 
of u, 5. q maps for several choices of 6, and 6. q maps for a sequence of temperatures. 



13 



1. (q, Lu) maps for a = 



Figure [7| shows intensity maps of Imx(q, u) for a sequence of hole densities 5 for the odd 
channel (left panels) and the even channel (right panels); the q scan direction is shown in 
the inset of Fig. IHl The threshold energy Uthiq) [Eq. (j221)] is also plotted with a gray dotted 
line. Figure El^a) shows a result in the vicinity of the AF instability. Strong intensity is seen 
well below ujt^{q), indicating collective particle- hole excitations with the upward dispersion 
described by Eq. fj25|) . This soft collective mode reflects the commensurate AF instability at 
6 ~ 0.064. The mode disperses as a function of q and mixes into the particle-hole continuum 
at sufficiently large distance from Q, where it becomes overdamped and thus faintly visible 
on the color scale in Fig. [7||^a). When 6 is slightly increased, the soft collective mode rapidly 
shifts to higher energies. It has a nearly fiat dispersion around q = Q for 5 = 0.08 [Fig.[7||^b)] 
and a clear downward dispersion for 6 = 0.12 [Fig. [7|^c)]. In the particle-hole continuum 
above ujtii{q), significant intensity is seen in Figs. [Z(b) and (c) with an upward dispersion. 
This can be interpreted as an overdamped collective mode near Wq^, as we shall discuss 
later. With further increasing 6, the collective mode is pushed up to higher energy and 
appears only close to u = UthiQ.) [Fig. Q^d)]. Above u^hiQ), the spectrum broadens and 
the structures become less clear. For 6 = 0.20 [Fig. El^e)], the collective mode does not 
appear and we just see that the continuum spectrum has strong intensity near q = Q and 
u = Uti^{Q), and a very broad featureless distribution above. As a function of 6, Figs. [7|^a)- 
(e) show that the spectral weight becomes larger at lower 6 (see values of the color map 
index) due to the proximity to the AF instability. 

The corresponding results for the even channel are shown in the right panels in Fig. [7| 
The overall features look the same as those of the odd channel. Close to the AF instability 
[Fig. El^f)], however, the collective mode retains the downward dispersion and does not show 
softening, refiecting the fact that the AF instability is driven by the odd channel. The 
collective mode is pushed up to higher energy with 6 [Figs. Wig) and (h)] and appears inside 
the continuum in Fig. [3^i) in the sense that Eq. ()2Hj) is still satisfied. The mode in the 
continuum is substantially overdamped but has the same |q — QP dispersion as Eq. (|^. 
The coefficient, however, is not given by Eq. ()26|) . which is valid only if Imxoi'i, ^^q^) = 0. 
Figure [7(j) looks similar to Fig. [7[i), but Eq. is not satisfied any longer and the strong 
intensity results from individual excitations. Compared to the odd channel, we see that 
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spectral weight of the even channel is overall suppressed and broadened. This difference 
comes from two effects: (i) the suppression of |t/(Q)| for the even channel near q = Q [see 
Eq. (j2(J|) ]. and (ii) different particle-hole scattering processes, namely interband scatterings 
for the odd channel and the intraband scatterings for the even channel. The former (latter) 
effect is dominant for low (high) 6. 

To see the intensity profile of Imx(q, u) more clearly, we replot Fig. [7|^c) by separating 
it into three energy regions in the left panel in Fig. |H1 with an optimized color map scale 
in each; the right panel is a map of |1 + J(q)Rexo(q, ^)\, which serves to quantify the 
collective character of magnetic excitations [see Eq. ()2H|l ]: the cross symbols indicate peak 
intensity positions of Imx(q, uj). For u < uJq^ = 0.403 J, the ingap collective mode has 
a downward dispersion. This collective mode is however realized only in a limited energy 
region. See the region where |1 + J(q)Rexo(q; ^) \ becomes zero in the right panel; at lower 
07,11 + J(q)Rexo(ci) ^) \ remains finite such that the peak in Imx(q, tu) does not correspond 
to a genuine collective mode. With further decreasing u the IC signals are substantially 
diminished, while the DIG signals continue down to zero energy. For uj > uj'q, we see another 
dispersive structure inside the continuum. This upward dispersion can be interpreted as an 
overdamped collective mode up to ~ 0.50 J in the sense that |1 + J(q)Rexo(q) ^)| becomes 
zero near the peak position. Note that this overdamped mode is not directly connected with 
the downward collective mode below uJq^. Although the overdamped mode seems to continue 
smoothly up to higher energy in the left panel, it changes into a peak of individual excitations 
as seen from the finite value of |1 + J(q)Rexo(q5 ^)\ in the right panel. 

In Fig. El we plot Imx(q, cu) at q = Q as a function of u for (a) 6 = 0.08 and 0.15, 
and (b) 6 = 0.20. The sharp peak for 6 = 0.08 corresponds to the collective mode; the 
continuum cannot be seen on the scale of Fig. Ofa). For 5 = 0.15 in the odd channel, we 
see both the collective mode at a; = 0.42J and the continuum above u = 0.45J; the two 
are not completely separated because of finite F and T. The peak at u = 0.49J for the 
even channel is a collective mode, but its peak width is much broader than that for the odd 
channel, indicating a much shorter life time of the mode. The odd channel spectrum for 
6 = 0.20 is similar to that for 6 = 0.15. The peak at = 0.38 J, however, is not a collective 
mode; Eq. ()23p is not satisfied any longer. Similarly the peak for the even channel results 
just from individual excitations. 

We summarize the peak position of Imx(q, tu) at q = Q, which we refer to as = uq, as 
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a function of S in Fig. [TUl ujq is equivalent to when Eq. is satisfied; the lower edge 
energy of particle-hole continuum, u;th(Q), is also plotted by the bold gray lines; we consider 
the region, 6 > 0.064, below which the AF instability takes place. With increasing doping, 
ujq increases rapidly for both channels. For the odd channel, ujq has a broad maximum 
around 5 = 0.15 and merges into the lower edge of particle- hole continuum; the particle- hole 
bound state is realized for 6 < 0.17 (solid circles in Fig. [TUl) . For the even channel the 
bound state is realized in the same doping region as for the odd channel (solid triangles), 
but appears inside the continuum in 0.12 < 5 < 0.17 accompanied by the discontinuous 
change of uJq^. This jump is due to the fact that the even channel consists of two different 
intraband scattering processes, within the bonding band and within the antibonding band. 

2. q maps for a = 

We next show q maps of Imx(q, u) around q = Q for a sequence of energies u for 

5 = 0.12. At low uj [Fig. fTlT a)]. strong spectral weight is localized in the region around 
q = (tt ± 2711], IT ± 2717]), which is due to the particle- hole excitations around the d-wave gap 
nodes. With increasing u, the scattering processes with q = (tt ± 27r?7, vr) and (tt, tt ± 2711]) 
begin to contribute and the strong intensity region forms a diamond shape [Fig. ITlT b)]. 
This diamond shrinks when uj increases further [Fig. ITTT c)]. starts being due to collective 
excitations [Fig. ITlT d)]. and is finally reduced to the collective commensurate peak at u; = 
cUq^ = 0.403 J [Fig.lTlTe)]. Above uJq^ [Fig.HH^f)], the DIG spectral weight becomes dominant 
with rather broad features. While the spectral weight difference between the DIG and IG 
positions is not sizable and the DIG structures are smeared out in a certain energy range 
[Fig. HHg)], the DIG signals seem to disperse outwards with increasing u [Figs. ITlTfVfh)] 
as expected from Fig. |S1 

The spectral weight distribution in q-space for other doping rates, S = 0.08 and 0.15, 
is shown in Fig. for uj < uJq^ (left panels) and uj > (right panels). For u < uJq^ 
we see that the strong weight forms a diamond shape independent of 6, and the IG signals 
become stronger than the DIG peaks for higher 6. For u > uJq^, on the other hand, the DIG 
structure is clearly seen at relatively low 6, while the spectral weight spreads broader with 

6 (see also Fig. [7j). 

To understand the intensity distribution in Fig. ^2 we write the imaginary part of Eq. (fTIHl 
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as 



Imx(q, ^) = n , 7/„at^„.. /„ . , r 7/„xt^... /„ . ' V^^) 



imxo(q, ^) 

[1 + J(q)Rexo(q, + [J(q)Imxo(q, uj) 
We have to consider both | J(q)Imxo(q, ^)\ and |1 + J(q)Rexo(q! ^)\, which we show in 

the left panel and the right panel in Fig. UHl respectively, at the same u as in Fig. [TTJ 

the structure of Imxo(q, ^) is essentially the same as that of J(q)Imxo(q7 Two typical 

energy scales, a;Q*^ and c<Jth(Q), characterize the spectral weight distribution of Imx(q, u). 

(1) For LU <t: uJo", we find |1 + J(q)Rexo(q, ^)\ > |-'^(q)ImXo(q7 ^)\ in Figs. CSta), 
(a'), (b), and (b'). Since the q dependence of Rexolq^ ^) is weak around q = Q we have 
Imx(q, u) oc Imxolq, that is, the weight distribution in Figs.[TTfa) and (b) is determined 
mainly by Imxo(q; 

(2) For u < cjq^ [Figs. \Wlc) and (c')], the minimum values of |1 + J(q)Rexo(q5 ^)|, 
become comparable with | J(q)Imxo(q; (x')l in a certain q region. From Eq. ()27|) . therefore, a 
strong weight of Imx(q uj) is determined by the minimum position of |1 + J(q)Rexo(q; ^)|; 
not by peak position of Imxo(q, Since Imxo(q, does not have an apparent structure 
in the region where |1 + J(q)Rexo(q; ^)\ takes minima, it depends on details which signal 
can be stronger, the IC or the DIG. In our model, the IC signals tend to be stronger than 
the DIG signals at higher 6. 

(3) At ~ Uq^, |1 + J(q)Rexo(q5 ^^)| can vanish for suitable q vectors [Figs. IT^ d') and 
(e')], leading to Imx(q, u) oc Imxo(q) uj)^^. Since Im^olq, ^) also can become very small 
in the same q region [Figs. [T^fd) and (e)], the strongest weight of Imx(q, a;) appears there, 
corresponding to a collective mode. 

(4) For Uq" ^ ^ ^^th{Q), |1 + <^(q)ReXo(q, ^^)| vanishes in the region where 
Imxo(q5 t^) has IC peaks [Figs. fT^ f) and (f)]. Since Imx(q, ou) oc lmxo{<i , when 
1 + J(q)Rexo(q, t^) = 0, the IC spectral weight is substantially diminished and Imx(q, u) 
shows Die peaks as seen in Fig. [TTff). 

(5) At u ^ c<Jth(Q), the IC peaks of | J(q)Imxo(q, merge into a commensurate peak as 
shown in Fig. fT^ g). | J(q)Imxo(q; ^^) \ does not have an apparent structure around the region 
where |1 + J(q)Rexo(q5 ^^)| vanishes [Fig.HSfg')]. Hence the strong weight of Imx(q,u;) has 
a similar distribution to that of 1 + J(q)Rexo(q, ^) ~ 0, and shows a "ring" in Fig. ITIT g). 

(6) For u > u;th(Q), Imxlq? ^) has strongest intensity in the region where |1 + 
J(q)Rexo(q; ^)\ ~ |-'^(q)IixiXo(q) ^^)| [compare Fig. [TTfh) with Figs. ITST h) and (h')]. 

The spectral distribution of Figs. fT^ a)-(d) corresponds to the cases (2), (4), (2), and (6), 
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respectively. 

3. T dependence for a = 

So far we have presented results for a fixed temperature T = O.OIJ. Now we discuss 
the T dependence of Imx(q, u). In Fig. [T31 we plot Imx(q, cu) at q = Q as a function 
of u) for several choices of T for 6 = 0.10 in the odd channel. The ci-wave gap disappears 
at Trvb = 0.139J for this doping. We see that the sharp peak of Imx(Q, u) survives 
as long as the d-wave gap exists. The peak energy is determined mainly by minimizing 
|l + t/(Q)Rexo(Q5 ^^)\- The T dependence of the peak energy is shown in the inset of Fig.lT^ 
it follows approximately the T dependence of the d-wave gap [see Fig. |2fa)]. With increasing 
T low-energy spectral weight proportional to uj becomes noticeable, and Imx(Q, uj) evolves 
smoothly toward the spectrum in the normal state (T > Trvb), where the sharp peak at 
finite cu disappears. 

The temperature dependence of q maps is shown in Fig. ^1 for several choices of T 
at 6 = 0.10 and u = 0.25J. Strong intensity of Imx(q, u) forms a diamond shape at 
low T. Although the diamond shape is retained with increasing T, the weight inside the 
diamond becomes stronger, leading to a dominant commensurate signal above T ~ 0.1 IJ. 
The commensurate signal thus dominates already at tempertures slightly below the onset 
temperature of (i-wave singlet pairing Trvb = 0.139 J. 

A (q, u) map of Imx(q, u) is shown in Fig. \TEi a) for T = 0.12J (< Trvb), where the 
energy region is separated into three with different color scales; the cross symbols mark 
highest weight positions along q = {qx, vr) and q)- The strongest weight appears at 

q = Q and uj ^ 0.26 J. Below this energy IC structures are not realized any longer and 
the spectrum shows a broad commensurate structure. However, in the high energy region 
{uj > 0.30 J) DIG signals still survive as well as IC signals. When T is further increased above 
^RVB [Fig. ll6r b)]. neither IC nor DIC structures are realized. The strong intensity is centered 
around q = Q and broadly spreads out. The right panels in Fig. are corresponding 
results for the even channel. Overall features look similar as for the odd channel except 
that the spectral weight becomes rather broadened. Although the peak positions form a 
rather complex pattern in Fig. ITHT c). an upward dispersion for 0.3 J ^ ^ 0.45 J can still 
be recognized together with a remnant of the downward dispersion for 0.2 J ^ uj < 0.25J. 
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4- uj dependence of q maps for a ^ 



In the following three subsections, we present results for Imx(q, cj) for a = —0.05 to 
discuss effects due to an in-plane anisotropy and their enhancement by (iFSD correlations. 

Figure El shows q maps of Imx(q, u) for a sequence of energies u for 6 = 0.12 and 
T = O.OIJ; the corresponding results for a = were shown in Fig. El At low u, the 
DIG signals appear, and form a high intensity region which is a bit distorted compared to 
Fig. ITlT a). due to the in-plane anisotropy. These DIG signals come from nodal scatterings. 
Note that the nodes of the singlet pairing deviate slightly from the diagonal direction because 
of 7^ Aj^ and finite A^ in Eqs. © and (jH)). Figures ITTT bVfd) are results for larger 
uj < uJq^. Appreciable weight appears on both the and qy axes and in this sense a nearly 
two-dimensional excitation spectrum is obtained even if the antibonding FS becomes open as 
in Fig. Analogous behavior was already found for the single layer t-J model in Ref. 
Effects of the anisotropy are seen in (i) the difference of the IG peak positions, namely 
(tt ± 27ii]x, tt) and (tt, tt ± 27criy) with rix 7^ rjy, leading to a deformed diamond shape, and 
(ii) the relative peak intensity difference of the two IG peaks. The latter however depends 
strongly on u and 6 even if a is fixed. For example, the IG peak at q = (tt, 7i±2TTriy) is a bit 
higher than that at q = (tt ± 2TTrix, tt) aX u = 0.20 J [Fig. ITTT b)]. while the order is reversed 
at uj = 0.35 J [Fig. fTTTc)]. Moreover, as we shall discuss in Sec.V.B., the peak intensity 
difference can depend strongly on details of the band structures and thus on models. On 
the other hand, the former is a robust feature of Imx(q, u). When the difference of the 
IG peak position, Arj = rj^ — rjy, is taken as a measure of the anisotropy of the IG peaks, 
we have Arj > for a < and in addition Arj oc a in good approximation at least up to 
\a\ ~ 0.1. Furthermore, the uj dependence of Arj is weak. In Fig. El we plot rjx and riy as a 
function of u. While 7]^ and riy depend on u, one finds At] > for all u < uJq^, and a weak 
u dependence of Arj except for the region of ~ uj'q where both rj^ and rjy tend to zero. 

Figure ITTT e) corresponds to the energy u = uJq^. While the strongest weight appears 
at q = Q for that energy, Imx(q, cUq*^) has an apparent elliptic distribution in the q plane. 
Figures E^f)-(h) are for u > uj'q. In Fig. ITTT f). the strongest weight appears along the 
diagonal direction, which is a remnant of the DIG signal for a = [Fig. . The strong 
weight distribution, however, forms a rectangular shape and shows a pronounced anisotropy. 
As energy is increased a little, Imx(q, uj) shows one-dimensional like signatures in Fig.Efg) 
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even with a small anisotropy (5%). To understand this one- dimensional pattern, we go back 
to Figs. Uni (g) and (g'). For a = —0.05, Imxolq, ^) still shows a commensurate signal but 
the spectral weight spreads out largely along the direction. Since |1 + J(q)Rexo(q; ^)\ 
vanishes around q = Q as shown in Fig. IT^ g'). the resulting Imx(q, uj) shows enhanced 
spectral weight along the qy direction, leading to Fig. ITTTg). With further increasing u, 
however, the spectral weight becomes much broader [Fig. fTTT li)] and the anisotropy becomes 
less apparent. For uj > Uq^ therefore the spectral weight distribution forms various patterns 
and is very sensitive to energy. This is in sharp contrast to the robust feature of the deformed 
diamond shape distribution for u < Uq", which enables us to define Ar] as a measure of the 
anisotropy of the magnetic excitations. 

5. 5 dependence of q maps for a ^ 

The left panels of Fig. El show q maps of Imx(q, u) for a sequence of doping concen- 
trations 6; the energy uj is below uj^q (actually below ujq for 6 = 0.20, where no collective 
mode exists, see Fig. ITUI). Appreciable intensity appears along a deformed diamond for all 
6. Although the renormalized band anisotropy decreases with increasing 6 [Fig. |2fb)], the 
anisotropy of Imx(q, u) is enhanced. To identify the influence of (iFSD correlations, we 
compare these results with the corresponding results for the bare anisotropy on the right 
panels of Fig. El The latter are calculated by switching off the rfFSD correlations, that 
is, we impose the same bare anisotropy a = —0.05, but use mean fields for a = 0. We 
see that the effect of the bare anisotropy on Imx(q, uj) is doping dependent even for fixed 
(doping-independent) a and becomes less pronounced at lower 6. In particular, the dia- 
mond is almost symmetric for 6 = 0.08 and 0.12 [Figs. fTW a') and (b')]. That is, the c/FSD 
correlations are particularly important to understand the anisotropy of Imx(q, uj) at lower 
S [Figs. imT a) and (b)]. This holds for both odd and even channels and similar results are 
obtained for the even channel. 

6. T dependence of q maps for a / 

Momentum space maps of Imx(q, u) for a sequence of temperatures are shown in the left 
panel of Fig. 1201 for 6 = 0.10 and uj = 0.25J. The deformed diamond shaped distribution is 
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rather robust against T [Figs- I^UT a) and (b)], although the spectral weight inside the diamond 
gradually increases with T, as we have seen in Fig. ^| At temperatures slightly below 
^RVB = 0.133 J, the commensurate signal becomes dominant [Fig. EUT c)] and a pronounced 
anisotropy appears in Fig. EUT d). When T increases further, the anisotropy is however 
reduced [Figs. EUT e) and (f)]. The temperature dependence of the anisotropy in Fig. QUI is 
not monotonically linked to the T-dependence of the band anisotropy. As seen in Fig. Efb) 
for 6 = 0.10, the band anisotropy at T = 0.13J is larger than that at T = 0.12J, but does 
not lead to a stronger anisotropy of Imx(q, u) at T = 0.13 J. 

To analyze the relevance of (iFSD correlations, we calculate Imx(q, u) also for isotropic 
mean fields and show the results in the right panel of Fig. 1201 We see that while (iFSD 
correlations contribute to the enhancement of the anisotropy of Imx(q, uj) at any T, such 
an effect becomes most pronounced at relatively high T as seen in Figs. EOfd) and (d'). 
While Fig. 1201 has been obtained for fixed u and 6, results for other parameter sets including 
the even channel show that the (iFSD correlations typically drive a pronounced anisotropy 
at relatively high T (~ Trvb)- 



V. DISCUSSION 



We have investigated Imx(q, uj) in the slave-boson mean-field scheme for the bilayer t-J 
model. We summarize our results through a comparison with inelastic neutron scattering 
data for YBCO. Experimental data are well summarized in Ref. 

S, Q- It should 



be kept in mind that in the slave-boson mean- field theory |43[ Trvb has to be interpreted 
as pseudogap crossover temperature T* in the underdoped regime, and as superconducting 
phase transition temperature in the overdoped regime of high-T^ cuprates. 



A. Prominent features of magnetic excitations in YBCO 



One of the most prominent features of mag netic excitations in YBCO is that IC peaks 



appear only below Tc or T' 
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121 Il3| which is in sharp contrast to LSCO where 
IC peaks were reported even at room temperature. This is captured by the present 

theory. As shown in Fig. the IC signals are obtained only at low T and change into 
a broad commensurate signal above T ~ Trvb- It was pointed out a few years ago that 
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the distinct behavior of LSCO is also captured by the present theoretical framework with 



different band parameters. 



45| 



Another central issue in YBCO is the so-cal 
a sharp magnetic signal at q = Q and finite uj. 



alternative scenarios for this resonance peak. 
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112l | the resonance is interpreted as a particle- hole bound state (Figs. El and [7j). While 



the resonance is particularly pronounced at low T Trvb), it starts to develop below 
^RVB (Fig. EH), that is, below T* in the underdoped regime and below in the oyerdo ped 
regime. This is consistent with experiments for optimally doped YBCO 



and YBa^CucjOfi 



for both odd 
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channels if we set J ~ 100meV.|79l| Although 



this agreement is based on tuning one of the band parameters, namely t/ J, it is remarkable 
that the 5 dependence of cuq*^ is reproduced within the present simple framework. We note 
that no genuine resonance was obtained for 5 > 0.18 in Fig. EH which however depends on 
the choice of F. 

As seen from Figs. [7| and |H1 the collective mode usually has a downward dispersion and is 
necessarily accornpanied by IC signals for uj < uj^q. This downward dispersion was actually 



observed 



Ly accoE 
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22j | The resonance at q = Q and the IC peaks observed in experiments 



are therefore understood as coming from the same origin. The resonance mode, however, is 
realized only in a limited energy region in the presence of finite damping and temperature, 
and loses its collective nature when it touches the continuum spectrum (Fig. Ej). 

The IC peaks are well defined along the cut with q = (g^., tt) or (tt, g^) for uj < cUq''. In 
the q plane, however, strong intensity weight appears on a diamond shaped region as we 
have seen in Figs. ITlT bVfd). Hence, the DIC peaks are also well defined along a cut with 
q = (g, q) / \pl. Although the possibility of DIC peaks was not discussed in the experimental 
literature, this diamond shaped distribution was actually observed. Q For uj <^ uJq^ the DIC 
signals become dominant [Fig. llir a)]. which is a robust feature comin g from th e rf-wave gap 
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113| but have not 



nodes. These DIC signals were already predicted about a decade ago, 
been detected in experiments. 

For UJ > lJq^, we have obtained an upward dispersion in Figs. [7| and |H1 especially at 
relatively low 6 and have interpreted it as an overdamped collective mode; note that this 
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mode is not directly connected with the downward collective mode in 00 < uJq^. An upward 
dispersion of DIG peaks for uj > uJq^ can be read off also from Figs. ^2 and ^| These 
DIG peaks are robust at least for ^ ^ ^^th(Q) as discussed in Figs. ITHT f) and (f), 
although they seem to continue up to higher uj than c^thlQ) in Figs. [7|and SI An upward 
DIG dispersion was recently reported in neutron scattering experiments. 0, 22, 

As seen in Figs. IT^ a) and (c), the upward dispersion in the high energy region is rather 
robust to temperature while in the low energy region temperature spoils the incommensurate 
structures in favor of a broad commensurate signal. This behavior was detected in the most 
.ecent.e„..o„scattenngexpe..ents.a 

The even channel behaves qualitatively similarly to the odd channel as seen from Fig. [7| 
The main differences are that its intensity is overall suppressed and it has a larger energy 



scale (characterized by u;q, see also Fig. HU)). These features o: 
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the even channel were actu- 



84 l85l | The experiments, 



ally reported in neutron scatting experiments. |2z 
however, indicate a much stronger suppression of the even channel than that obtained in our 

n 

calculation, especially in the low uj region. In fact, Fong et a/.|S^ concluded that the even 
channel was fully gapped even at T = 200K. Although the even channel could be suppressed 
more by taking a larger value of in Eq. (j^UI) . such a "gap" feature is not captured by the 
present theory. 

While the present theory captures prominent features of magnetic excitations in YBCO, 
the line shapes of Imx(Q, uj) in Fig. 121 are different from experimental data. Q,S,Q The 
experimentally observed peak structure, which is interpreted as resonance, is much broader. 
From the experimental viewpoint it is not evident that the resonance is really situated below 
the continuum, although we have interpreted it as an ingap collective mode especially for 
the odd channel (Fig. ITUI) . Several experimentalists jl^ lioi l84 | extract the gap from the 
low-energy tail of the resonance, such that the gap becomes smaller than the resonance 
energy by definition. From the present theoretical viewpoint, however, the gap should be 
defined from the lower edge of the continuum. 



B. Anisotropy of magnetic excitations 



Hinkov et a/. [54] presented a clear q map of magnetic spectral intensity for untwinned 
YBa2Gu306.85, which showed that the peak intensity at q = (7r±27rr7^, vr) is larger than that 
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at q = (vr, vr ± 2Tiiriy) . This map was obtained by fitting the exp erime ntal data to putative 



(5-peaks in q space smeared by a spectrometer resolution function. |114| In the present model 
(see Fig. IT^ . the peak intensity difference for uj < uo^^ is less sizable below 6 ~ 0.15 
than the difference extracted from the fit to the experimental data. Hinkov et a/.js^ also 
reported that the q-integrated spectral weight difference between along q = (tt ± 2Txrix., tt) 
and q = (vr, vr ± 2TT'r]y) became apparently larger with decreasing u below the resonance 
energy. This feature is not manifest in our results shown in Figs. ITTT cVfe). 

These discrepancies may result from details of the band structure. LDA calculations 
suggest 0] that seemingly small effects, such as long-range hopping integrals beyond third 
neighbors and the buckling of oxygen atoms, may affect the dispersion near the Fermi level 
without necessarily leading to strong shifts of the Fermi surface. The recently observed 
kink structure of the band dispersion near the Fermi level 115l | should also be taken into 



account. The influence of band structure features on the IC peak intensity difference follows 

n 

also from theoretical studies. Zhou and Li|58] showed that the direct coupling to the chain 
band contributes strongly to the peak intensity difference. Two recent phenomenological 
calculations with different band parameters also revealed the sensitivity to band structure 
details: in one case 22% band anisotrop yiSftl was required to account for the peak inten- 
sity difference reported by Hinkov et a/.jS^, in the other only 6%. There may also be 
ambiguities in the analysis of the experimental data, since the anisotropy of the IC peaks 
was deduced after a background correction, and the larger anisotropy was reported when 
the background became dominant, that is when the magnetic signals became very weak. If 
there is a magnetic signal hidden in the background, one might overestimate the IC peak 
intensity difference. 

Since the anisotropy of the IC peak intensity is sensitive to details of the model and 
to ambiguities in the experimental data analysis, it seems advantageous to consider the 
anisotropy of the incommensurability. At] = 77^ — r^^, as a more robust measure of the 
anisotropy of Imx(q, uj) at low T for uj < uj'q. We have consistently obtained Ar/ ^ for 
a ^ for all studied parameter sets, and moreover At] turns out to be proportional to a in 
good approximation at least up to |a;| ~ 0.1 without appreciable uj dependence except near 
cu^*^ (Fig. CHI). Our resuh, Ar] ^ 0.02 - 0.03 for a = -0.05 (Fig. CHI), is comparable with 
experim ental data for YBa2Cu306.6 and not inconsistent with the data for YBa2Cu306.85 
also. 
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The negative sign of a corresponds to < ty, that is, the opposite of what one 
would expect from the in-plane lattice constant difference. In agreement with LDA band 
calculations, IstI this implies that the chains are crucial to the in-plane anisotropy in YBCO 
and thus to the understanding of the anisotropy of Imx(q, uj). The importance of chain ef- 



fects was already discussed by Zhou and Li 



in a different context. They analyzed direct 



coupling between CuO chains and Cu02 planes. Such a direct coupling, however, should be 



small enough to ensure that Imx(q, u) is wel 



sistent with the observed modulation 



.0 



13 



characterized by a bilayer model to be con 
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17 



19 



79 



80 



We have therefore 



interpreted chain effects as mainly renormalizing the in-plane band anisotropy. We note that 
an LDA calculation |j82i] predicts that chain effects are not important to the anisotropy of t 
in the double chain compound YBa2Cu408, for which we thus expect > ty. Our theory 
then predicts Arj = rj^ — rjy < for uj < uJq^ at low T. 

For u > u'q we have found a pronounced spectral weight anisotropy, which however 
strongly depends on u [Figs. IT7r f)-(h)]. Stock et al^^ performed high energy neutron 
scattering experiments for partially untwinned YBa2Cu306.5 and obtained q maps of the 
spectral weight. They reported nearly isotropic spectral weight distribution. However, their 
data were integrated over intervals ±7.5meV along the energy axis, which may smear out the 
predicted anisotropy in Figs. IT7r f)-fh). In addition, anisotropics may be underestimated if 
the untwinning of the crystal is not complete. On the other hand, Hinkov et a/.js^ observed 
an anisotropy of IC peaks at several choices of energy for uj > ujq^ in almost fully untwinned 
YBa2Cu306.6, while full q maps have not yet been obtained. Further neutron scattering 
data for u > would be useful. 

At relatively high T we have found an enhanced anisotropy of Imx(q, u) [Figs. OUT c) 
and (d)]. This anisotropy is characterized by the difference of the broad commensurate 
peak width between the direction and the qy direction. This anisotropic peak width 
was actually observed in the most recent experiment for YBa2Cu306.6 in the pseudogap 
regime. |3] 

C. Relevance of tZFSD correlations 



A sizable anisotropy of Im^f 
YBa2Cu306.5Q and YBa2Cu306 



''6.6- 



uj) was reported for underdoped materials, 
\5^ To understand this behavior, the bare band 
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anisotropy, without an enhancement due to c/FSD correlations, is not sufficient. The mag- 
netic anisotropy due to the bare band structure anisotropy decreases at lower 6 for fixed a, 
as seen in the right panels of Fig. ^1 One cannot expect a larger a at lower 6 since the 
crystal structure changes from orthorhombic to tetragonal for y < 6.4, and is accompanied 



by the disappearance of the CuO chains. 



90| On the other hand, the presence of dPSD cor- 



relations provides a natural explanation for the observed anisotropy at low 6; the correlation 
effects can yield an enhanced anisotropy of magnetic excitations at low 6 fFig. I19|). While 
we have treated (iPSD correlations in the slave-boson mean-field approximation to the t-J 
model, they were shown to be present also in a recent exact diagonalization study, jsli The 
anisotropy of Imx(q, uj) in optimally doped YBa2Cu306.85|3; other hand, can be 

understood qualitatively already from the bare band anisotropy [Figs. lTW c'l and (d')], which 
is however further enhanced by dPSD correlations as seen in Figs. ITW c) and (d). 

(iFSD correlations also provide a natural scenario to account for the observed enhanced 
anisotropy of Imx(q, uj) in the pseudogap phase of underdoped YBa2Cu3O6.6-0| As seen 
in Figs. EDI (c), (c'), (d), and (d'), the anisotropy of Imx(q, u) is substantially enhanced 
compared with the bare anisotropy effect especially at relatively high temperatures. We note 
that the obtained anisotropy at high T in the left panel of Fig. 1201 may be underestimated 
in the present calculation since we have assumed an isotropic hopping amplitude of bosons 
(&J&r') = ^- Above the superconducting transition temperature, the bosons are not really 
condensed at the bottom of the band and (&J&r') can become anisotropic, contributing thus 
to an enhancement of the anisotropy of magnetic excitations. 

Due to dFSD correlations, we can expect a relatively large anisotropy at low 6 and high T 
in orthorhomb ic YB CO. A large anisotropy was actually reported in resistivity measurements 
a few yeas ago. 116^ The data, however, were interpreted in terms of partial spin-charge stripe 
order, often referred to as (electronic) nematic order. 3| Although the nematic order has the 
same symmetry as the (iFSD, namely orientational symmetry breaking of the square lattice. 



the underlying mechanism is different. ciFSD correlations come 



rom forward scattering 



processes of quasiparticles close to the FS near (vr, 0) and (0, 7r),|44,Q while a spin-charge 
stripe requires correlations with a large momentum transfer such as ant if erroma gnet i sm. 
Our scenario offers a different route to understand the aniso tropy o f the resistivityj'll6l] as 



) tropy o i tne 

10Q Q 



well as the anisotropy of the magnetic excitation spectrum |13l. |5J, |55[ on the basis of a 
microscopic model calculation. 
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D. Fermi surface and superconducting gap anisotropy 



Because of dFSD correlations, the FS of the antibonding band can easily open in presence 
of a small bare band anisotropy as shown in Fig. EJ The Fermi surface topology depends 
in part icular on the values of a and tzjt. Further efforts to determine the FS shape in 



YBCO 
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118| will serve to extract these parameters, which then gives insights on chain 
band effects as well as on the present scenario for the anisotropy of magnetic excitations. As 
shown in Fig. the hole density of the antibonding band is significantly larger than that 
of the the bonding band. This self-consistent result may be useful for mapping the FS in 
future ARPES studies. 

Although the inner FS can change its topology with a, such a topology change does not 
strongly affect the anisotropy of Imx(q, lS). We have not observed an enhanced Ar] when 
the inner FS opens, but a simple linear relation Ar] oc a in good approximation at low T for 

uj < cJ^^ 

The present self-consistent calculation yields A^; > |Aj,| for a < 0, that is, Ag > 
[Fig. |2fa)], which partially compensates the bare anisotropy of Jrr' [see Eq. (0)]; note that 
A^ is the RVB pairing amplitude and the true gap magnitude is given by J^A^. Hence 
we obtain Jx^x < {Jy^yl with an anisotropy a bit smaller than \2a\, about 7% for 6 = 
0.10 — 0.20 and a = —0.05. The relation Jx^x < \Jy^y\ "was actually r epor ted in a Raman 
scattering experiment for the optimally doped and overdoped YBC0.(ll9| Note that this 
gap anisotropy does not lead to a gap anisotropy between k = (tt, 0) and (0, vr), which 
comes from A^ in the present theory, as seen from Eqs. (0) and ©. While the magnitude 
of A^ might look sizable in Fig. |21(c) especially for high 6, A^ is multiplied by in Eq. (jUj). 
The .e.ult.ng a,n.ot.opy between A,.,.,^.d A,.,„,„ becomes very =n,all, le. than 1% 
for our parameters. Recent ARPES datajll8l| reported an anisotropy of about 50%, which 
cannot be understood in the present model. 



VI. CONCLUSION 



We have performed a comprehensive analysis of the bilayer t-J model in the slave-boson 
mean-field scheme. After determining the mean fields self-consistently, we have calculated 
the dynamical magnetic susceptibility in a renormalized RPA. We have shown that promi- 
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nent features of magnetic excitations of YBCO are captured by this scheme: (i) the IC and 
Die signals only for T < Trvb, (ii) the collective mode for T < Trvb, (iii) its downward 
dispersion for uj < cUq*^, (iv) the 5 dependence of cJq'^, (v) the overdamped collective mode 
for Lo > Wq*^ with an upward dispersion, (vi) robustness of this high energy dispersive feature 
to T, and (vii) spectral weight suppression of the even channel. 

In particular, the present theory, which includes c/FSD correlation effects in a self- 
consistent manner, provides a natural scenario to understand the observed anisotropy of 
Imx(q, ijj): (i) appreciable anisotropy also for low 5, and (ii) enhanced anisotropy at rela- 
tively high T (~ Trvb)- 

Although magnetic excitations of high-Tc cuprates are often discussed in terms of spin- 
charge stripes after the proposal by Tranquada et a/.,|4S] the present comprehensive study 
indicates that conventional particle-hole scattering processes around the FS are essential 
to magnetic excitations. An additional insight from this study is the importance of dPSD 
correlations for the anisotropy of Imx(q, uS). c/FSD correlations are due to forward scattering 
interactions, which were so far ignored in most theories. 
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FIG. 1: Notation for superexchange couplings and transfer integrals in the bilayer square lattice 
model. 
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FIG. 2: Mean fields xo, ^d, and Xz for a = 0: (a) T dependence for 6 = 0.10, and (b) 6 dependence 
at T = O.OIJ. 
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FIG. 3: (a) T dependence of Xd, ^s, and for 6 = 0.10 and a = —0.05; and A^ are multiplied 
by 4; note that Xd < because of a < 0. (b) T dependence of anisotropy of the renormalized band, 
2At/Aio, for several choices of 5. (c) 5 dependence of A^ and A^ at T = O.OIJ. 
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FIG. 5: Hole density in the antibonding band {6a) and the bonding band {6^) as a function of total 
hole density 6; note that 6 = {6a + 6h)/2. The inset shows the ratio 6a/ 6. 
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(^x,Jt) (jt,Jt) {q,q)/i2 



FIG. 6: Schematic plot of magnetic excitation spectra Imx(q, w); the q directions are shown in the 
inset. The hatched region represents the continuum spectrum and its lower edge is a;th(q)- The 
bold solid line oj^^^ represents a typical dispersion of the collective mode. The mode is softened 
near the AF instability leading to an upward dispersion as shown by the dotted line. 
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(0.4II,II) (O.Slt,O.Slt) (0.4jt,u) (O.Sji.O.Si) 

FIG. 7: (color) (q, uj) maps of Imx(q, w) for a sequence of doping concentrations 5 at T = O.OIJ 
and a = for both odd (left panels) and even (right panels) channels; the gray dotted line is the 
threshold energy a;th(q); the q scan directions are shown in the inset of Fig. El {OAtt, vr) < q < 
(vr, vr) and (vr, vr) > q > (0.5vr, 0.5vr). 
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FIG. 8: (color) (q, lo) maps of Imx(q, cj) (a) and |1 + J(q)Rexo(q, w)| (b) for 6 = 0.12, T = O.OIJ, 
and a = in the odd channel; q-scans as in Fig. [7| The cross symbols represent the highest weight 
positions of Imx(q, o->) along q = {q^, vr) and -j^{q, q), respectively. 




FIG. 9: uj dependence of Imx(Q, uj) for S = 0.08 and 0.15 (a), and 6 = 0.20 (b) for the odd and 
even channels; the results for 6 = 0.08 are multiplied by 0.2. 
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FIG. 10: Doping dependence of wq for the odd and even channels at T = O.OIJ and a = 0; the 
superscript "res" indicates that Eq. is satisfied at q = Q and uj = ljq; u;th(Q) is also plotted 
with bold gray lines for both channels. 
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FIG. 11: (color online) q maps of Imx(q, uj) for a sequence of energies ui in 0.67r < gx,^?/ ^ l-47r 
for 6 = 0.12, r = O.OIJ, and a = in the odd channel. 
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FIG. 12: (color online) q maps of Imx(q, uj) in O.Gtt < qx,qy < l-47r for 6 = 0.08 (upper panels) 
and 0.15 (lower panels) for T = O.OIJ, and a = in the odd channel. The energy is set below 
(above) iVq^ in left (right) panels: uj'q = 0.245J and 0.421J for 6 = 0.08 and 0.15, respectively. 
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FIG. 13: (color online) q maps of | J(q)Imxo(q; ^)\ (left panel) and |1 + J(q)Rexo(q> ^)\ (right 
panel) for a sequence of energies uj in O.Gvr < qx,Qy < IAtt for 5 = 0.12, T = O.OIJ, and a = in 
the odd channel. In (c')-(g'), the maps are restricted to q with |1 + J(q)Rexo(q, '^)| < 0.5 to get 
a better contrast in the interesting region. 
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FIG. 14: (color online) oj dependence of Imx(Q, oo) for several choices of T for 5 = 0.10, and a = 
in the odd channel; here Trvb = 0.139J; for T = 0.05J and O.IOJ, Imx(Qi ^) is multiplied by 0.3 
and 0.5, respectively. The inset shows the T dependence of the peak energy of Imx(Q) 
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FIG. 15: (color online) q maps of Imx(q, uj) for a sequence of temperatures in O.Gtt < Qx, Qy < 1.47r 
for 6 = 0.10, Lo = 0.25 J, and a = in the odd channel; here Trvb = 0.139 J. 
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FIG. 16: (color) (q, u) maps of Imx(q, cj) at T < Trvb = 0.139J (a) and (c), and T > Trvb (b) 
and (d) for 6 = 0.10 and a = 0; the left (right) panels are for the odd (even) channel; q-scans as in 
Fig. [7| The cross symbols in (a) and (c) represent the highest weight positions along q = {q^, vr) 
and -j=iq, q). 
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FIG. 17: (color online) q maps of Imx(ci; ^) for a sequence of uj in O.Gtt < Qx-^Qy < l-47r for 6 = 0.12, 
T = O.OIJ, and a = —0.05 in the odd channel. 
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FIG. 18: Energy dependence of r]x and rjy for 6 = 0.12, T = O.OIJ, and a = —0.05 in the odd 
channel. The jump of r]x and r/y at w ~ 0.3 J is due to fine peak structures of Imx(q; 
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FIG. 19: (color online) q maps of Imx(q, for a sequence of doping concentrations 5 for T = 
O.OIJ, and a = —0.05 in the odd channel; lo is chosen below cOq^; the right panels show the results 
obtained from the bare anisotropy without the enhancement due to dFSD correlations; q is scanned 
in O.Gvr < qx,Qy < lAir except for the panels for 5 = 0.20 where O.Svr < qx,Qy < l-5vr is taken. 
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FIG. 20: (color online) Left panel: q maps of Imx(q7 ^) for a sequence of T in O.Gvr < q^^ Qy < 1.47r 
for 5 = 0.10, w = 0.25 J, and a = —0.05 in the odd channel; here Trvb = 0.133 J. Right panel: 
corresponding results for the bare anisotropy, that is, without dFSD correlations; Trvb = 0.139J 
in that case. 
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